MPI-Ph/96-57 
TUM-T31-94/96 
HD-THEP-96-23 
hep-ph/9608365 

<6 
o\ . 

A 1996 Analysis of the CP Violating Ratio e'/e^ 

wo: y & / 

Andrzej J. Buras 1 ' 2 , Matthias Jamin 3 , Markus E. Lautenbacher 1 

ON 

1 Physik Department, Technische Universitat Miinchen, 
D-85748 Garching, Germany 
■ 2 Max-Planck-Institut fur Physik - Werner-Heisenberg-Institut, 

m ■ 

qq ' Fdhringer Ring 6, D-80805 Miinchen, Germany 

\ 3 Institut fiir Theoretische Physik, Universitat Heidelberg, 

^3 
a: 

> 

• r- 1 

X 

We update our 1993 analysis of the CP violating ratio e'/e in view of 
the changes in several input parameters, in particular the improved value 
of the top quark mass. We also investigate the strange quark mass m s 
dependence of e'/e in view of rather low values found in the most recent 
lattice calculations. A simple scanning of the input parameters within 
one standard deviation gives the ranges: — 1.2 • 10 -4 < e'/e < 16.0 • 10 
and < e'/e < 43.0 • 10" 4 for m s (m c ) = 150 ± 20 MeV and m s (m c ) = 
100 ± 20 MeV respectively. If the experimentally measured numbers and 
the theoretical input parameters are used with Gaussian errors, we find 
e'/e = (3.6 ± 3.4) • 10" 4 and e'/e = (10.4 ± 8.3) • 10~ 4 respectively. We 
also give results for ImV t *Vtd- Analyzing the dependence of e'/e on various 
parameters we find that only for m s (m c ) < 100 MeV and a conspiration 
of other input parameters, values for e'/e as large as (2 — 4) • 10~ 3 and 
consistent with the NA31 result can be obtained. 
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The measurement of e'/e at the 10 -4 level remains as one of the important 
targets of contemporary particle physics. A non-vanishing value of this ratio 
would give the first signal for the direct CP violation ruling out the superweak 
models. The experimental situation on Re(e'/e) is unclear at present. While 
the result of NA31 collaboration at CERN with He(e'/e) = (23 ± 7) ■ 10" 4 
clearly indicates direct CP violation, the value of E731 at Fermilab, Re(e'/e) = 
(7.4 ± 5.9) ■ 1CT 4 0, is compatible with superweak theories |§ in which e'/e = 0. 
Hopefully, in about two years the experimental situation concerning e'/e will be 
clarified through the improved measurements by the two collaborations at the 
1CT 4 level and by the KLOE experiment at DA$NE. 

There is no question about that the direct CP violation is present in the 
standard model. Yet, accidentally it could turn out that it will be difficult to see 
it in K — > Tin decays. Indeed in the standard model e'/e is governed by QCD 
penguins and electroweak (EW) penguins. In spite of being suppressed by a/a s 
relative to QCD penguin contributions, the electroweak penguin contributions 
have to be included because of the additional enhancement factor ReAo/Re^ = 
22 relative to QCD penguins. With increasing m t the EW penguins become 
increasingly important 0,|t|, and entering e'/e with the opposite sign to QCD 
penguins suppress this ratio for large m t . For m t ~ 200 GeV the ratio can even be 
zero ||. Because of this strong cancelation between two dominant contributions 
and due to uncertainties related to hadronic matrix elements of the relevant local 
operators, a precise prediction of e'/e is not possible at present. 

In spite of all these difficulties, a considerable progress has been made in this 
decade to calculate e'/e. First of all the complete next-to-leading order (NLO) 
effective hamiltonians for AS = 1 AS = 2 and AB = 2 | are 

now available so that a complete NLO analysis of e'/ e including constraints from 
the observed indirect CP violation (ex) and the B® — B® mixing (Am^) is 
possible. The improved determination of the V u b and V c b elements of the CKM 
matrix |12|,^3| , and in particular the determination of the top quark mass m t [[14 



had of course also an important impact on e'/e. The main remaining theoretical 
uncertainties in this ratio are then the poorly known hadronic matrix elements of 
the relevant QCD penguin and electroweak penguin operators, the values of the 
non-perturbative parameters Bk and ^B B F B and as stressed in pj the values of 



m s and Ajj-g. 



In 1993 we have presented a detailed NLO analysis of e'/e |7| using all in- 
formation available at that time. A similar NLO analysis has been made by the 



Rome group fi~5 [. In 1995 the latter group has updated their analysis to predict 



a very small ratio e'/e = (3.1 ± 2.5) ■ 10~ 4 essentially consistent with the 
superweak scenario. An analysis of e'/e with different treatments of hadronic 
matrix elements can be found in [17-19] and will be briefly discussed below. 

The purpose of the present letter is to update our 1993 analysis and to confront 
our new result with the one of the Rome group [16|] and with fll7}-|l9fl. 

Let us list the main new ingredients of our present analysis compared to the 
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previous one: 



The value of the top quark mass from CDF and DO ITI] . 



Updated values for the elements of the CKM matrix such as V u b and V c b 

mm 



m B ) coming from most recent lattice 
calculations, 



New values of the strange quark mass 



and QCD sum rule [23-2 



Improved value for a s (M z ) p6| , for which we take a s (M z ) = 0.118 ± 0.005 



corresponding to A- 



MS 



325 ± 80 MeV, 



• The inclusion of NLO corrections to the QCD factor 773 |Tl| entering the top- 
charm contribution in the AS* = 2 effective hamiltonian, the last previously 
missing ingredient of a complete NLO analysis of e'/e, and 

• Two distinct analyses of theoretical and experimental uncertainties. 

Let us first recall the basic formulae used in our new analysis, referring fre- 
quently to our 1993 paper []7J where further details can be found. In we have 
analyzed e'/e in the Standard Model including leading and next-to-leading log- 
arithmic contributions to the Wilson coefficient functions of the relevant local 
operators Imposing the constraints from the CP conserving K — > tttt data 

on the hadronic matrix elements of these operators we have given numerical re- 
sults for e'/e as a function of Aj^g, m t , m s and two non-perturbative parameters 
B^ 2 ^ and B^ 2 ^ which cannot be fixed by the CP conserving data at present. 
These two parameters are defined by 



(Q 6 (m c )) = 4 1/2) (Q 6 (m c )} ( o 



(vac) 



(Q 8 (m c )) 2 = Bf /2) (Q 8 (m c ))t c) 



where 
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are the dominant QCD and electroweak penguin operators, respectively. The 
subscripts on the hadronic matrix elements denote the isospin of the final tttt- 
state. The label "vac" stands for the vacuum insertion estimate of the hadronic 



matrix elements in question for which B$ 



(1/2) 



B 



(3/2) 



1. The same result is 



found in the large N limit [p?| ,|28 



B 



(1/2) 



1.0 ±0.2 29,30 and B. 



1(3/2) 



Also lattice calculations give similar results 
= 1.0±0.2 [HHH, Bf /2) = 0.81(1) 0. We 



have demonstrated in that in QCD the parameters B^ 2 ^ and B^ 1 *' depend 
only very weakly on the renormalization scale /1 when /i > 1 GeV is considered. 
The /1 dependence for the matrix elements in (|l|) is then given to an excellent 
accuracy by l/m^(fi) with m s (fi) denoting the running strange quark mass. The 



,(3/2) 
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scale fi = m c in ([I]) is a convenient choice for the extraction of matrix elements 
from the CP conserving data. 

At this point it seems appropriate to summarize the present status of the value 
of the strange quark mass. The most recent results of QCD sum rule (QCDSR) 
calculations p3|-p5| obtained at \i — 1 GeV correspond to m s {m c ) = 170 ± 
20 MeV with m c = 1.3 GeV. The lattice calculation of finds m s (2 GeV) = 
128 d= 18 MeV which corresponds to m s (m c ) = 150 ± 20 MeV, in rather good 
agreement with the QCDSR result. This summer a new lattice result has been 
presented by Gupta and Bhattacharya ||21|| . In the quenched approximation they 
find m s (2 GeV) = 90 ± 20 MeV corresponding to m s (m c ) = 105 ± 20 MeV. 
For nf = 2 the value is found to be even lower: m s (2 GeV) = 70 ± 15 MeV 
corresponding to m s (m c ) = 82 ± 17 MeV. Similar results are expected to come 
soon from the lattice group at FNAL [2^]. Certainly these results are on the low 
side of all strange quark mass determinations. 

Now, as an average from quenched lattice QCD, we can take m s (m c ) = 130 ± 
20 MeV and when averaged with the QCD sum rule value, we arrive at m s (m c ) = 
150 ± 20 MeV. This will be one of the ranges for the strange quark mass to be 
used in our analysis. On the other hand we cannot exclude at present that 
the ultimate values for m s will be as low as found in the most recent lattice 
calculations |21,^2 . In order to cover this possibility we will also present results 



for m s (m c ) = 100 ± 20 MeV. In table [I] we provide the dictionary between the 
values of m s normalized at different scales. To this end we have used the standard 
renormalization group formula at two-loop level with = 325 MeV. 



Table 1: The dictionary between the values of m s in units of MeV normalized at 
different scales with m c = 1.3 GeV. 



m s (m c ) 


75 


100 


125 


150 


175 


m s (2 GeV) 


65 


86 


108 


129 


151 


m s (l GeV) 


87 


116 


144 


173 


202 



It should also be remarked that the decomposition of the relevant hadronic 
matrix elements of penguin operators into a product of Bi factors times 1/mj. 
although useful in the 1/N approach is unnecessary in a brute force method like 
the lattice approach. It is to be expected that the future lattice calculations 
will directly give the relevant hadronic matrix elements and the issue of m s in 
connection with s'/e will effectively disappear. 

The details of the calculation of the Wilson coefficient functions as well as the 
determination of the hadronic matrix elements from the CP-conserving data can 
be found in and will not be repeated here. We will rather present an update of 
the analytic formula for e'/ e of ref . |34j which to a very good accuracy represents 
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our numerical analysis. This analytic formula exhibits the dependence of e' / e on 
m t , m s , Aj^L 5g 1//2 ' ) and B^ 2 \ It is given as follows: 



e'/e = lm\ t F(x t ) 



(3) 



where 



F(x t ) = P + Px X {x t ) + P Y Y {x t ) + P z Z {x t ) + P E E {x t ) 



(4) 



and 



ImA* = lmV t *V td = \V uh \ \V ch \ sin 5 = r] A 5 A 2 



(5) 



in the standard parameterization of the CKM matrix ]35| and in the Wolfenstein 
parameterization fl36"f , respectively. 

The m t -dependent functions in (fl) are given by 



X {x t ) 
Y (x t ) 
Zo(x t ) 



x t 



Xt + 2 3xt 
x t -l 
x t -4 

+ 
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(x t - l) 2 
2>x t 



x t -l (x t - l) 2 



In x t 
In Xt 



1 , 18x 4 - 163a;? + 259xf - 108x t 
~z^x t + l - 7T ^ ^ - + 

y 



144 (x t - 1 



32x 4 - 38x 3 - 15x 2 + 18x t 



72{x t - l) 4 



In x t 



Eo(x t ) 



z?(15 - 16x t + 4x?) , 

In x t + ; — ^ In x t + 

6(1 -x t ) 4 



x t (18 - llxt - x 2 ) 
12(1 -x t ) 3 



(6) 



with Xt = m 2 /M^. In the range 150 GeV < m t < 200 GeV these functions can 
be approximated to better than 1% accuracy by the following expressions 



X (x t ) = 0.660 a^ 575 
ZJxt) = 0.175 x 



0.93 



y ( z< ) = 0.315 x° t 78 
E (x t ) = 0.564 x t -°- 51 



(7) 



The coefficients Pj are given in terms of P, 
and mJm r ) as follows 



(1/2) 
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P, 



(1/2), 



m c ), P 8 



(3/2) 



P. 



(3/2), 



P 



(o) , 



158 MeV 



m s (m c ) + m d (m c ) 



(8) 



The Pi are renormalization scale and scheme independent. They depend however 
on A™. In table we give the numerical values of r 



(0) „(6) 



and r) for different 



values of A£L at u = m c in the NDR renormalization scheme. The coefficients r 



(0) 



,(6) 



MS 



and r| depend only very weakly on m s (m c ) as the dominant m s dependence 



has been factored out. The numbers given in table |2] correspond to m s (m c 



150 MeV. However, even for m s (m c ) ph 100 MeV the analytic expressions given 
here reproduce our numerical calculations of e'/e to better than 4%. For different 
scales fj, the numerical values in the tables change without modifying the values 
of the P^s as it should be. To this end also Pg and B^ 2 ^ have to be modified 
as they depend albeit weakly on \i. 

Concerning the scheme dependence only the r coefficients are scheme de- 
pendent at the NLO level. Their values in the HV scheme are given in the last 
row of table |2|. The coefficients r i; i = X,Y, Z, E are on the other hand scheme 
independent at NLO. This is related to the fact that the m t dependence in e'/e 
enters first at the NLO level and consequently all coefficients T{ in front of the 
m t dependent functions must be scheme independent. 

Consequently, when changing the renormalization scheme one is only obliged 
to change appropriately Pg and B^ 2 ^ in the formula for P m order to obtain 
a scheme independence of e'/e. In calculating Pj where i ^ 0, B^ 2 ^ and B^ 2 ^ 
can in fact remain unchanged, because their variation in this part corresponds to 
higher order contributions to e'/e which would have to be taken into account in 
the next order of perturbation theory. 

For similar reasons the NLO analysis of e'/e is still insensitive to the precise 
definition of m%. In view of the fact that the NLO calculations needed to extract 
ImAt (see below) have been done with m t = mj(mt) we will also use this definition 
in calculating F(xt) ■ The value for m t (m t ) corresponding to the average pole mass 
m poi e = 175 ± 6 GeV from CDF and DO 0, is rn t {m t ) = 167 ± 6 GeV. In what 

follows m t stands always for m t (mt). 



Table 2: AS = 1 PBE coefficients for various A^g in the NDR scheme. The last 
row gives the r$ coefficients in the HV scheme. 





A^ = 245MeV 


A^ = 325MeV 


Ag^ = 405 MeV 


i 


r (°) 


r (6) 


r (8) 
' i 


r (o) 


r («) 


r (8) 
' i 


r (0) 


r (6) 


' i 





-2.674 


6.537 


1.111 


-2.747 


8.043 


0.933 


-2.814 


9.929 


0.710 


X 


0.541 


0.011 





0.517 


0.015 





0.498 


0.019 





Y 


0.408 


0.049 





0.383 


0.058 





0.361 


0.068 





Z 


0.178 


-0.009 


-6.468 


0.244 


-0.011 


-7.402 


0.320 


-0.013 


-8.525 


E 


0.197 


-0.790 


0.278 


0.176 


-0.917 


0.335 


0.154 


-1.063 


0.402 





-2.658 


5.818 


0.839 


-2.729 


6.998 


0.639 


-2.795 


8.415 


0.398 



The inspection of the table shows that the terms involving rg and r z 
dominate the ratio e'/e. The function Zo(xt) representing a gauge invariant 
combination of Z°- and 7-penguins grows rapidly with m t and due to r% < 
these contributions suppress e'/e strongly for large m t 0,0. 
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In order to complete our analysis of e' / e we need the value of lm.\ t . To this end 
we will use the standard expression for Ek describing the indirect CP violation in 
K —>■ 7T7T and the corresponding expression for Am Bd which describes the B d -B d 
mixing. Since these expressions are by now well known we will not repeat them 
here. They can be found in ||37|| . We list here only the QCD factors rji (i = 1, 2, 3) 
and t}b relevant for e k and x d respectively: 



m = 1.38 7]2 = 0.57 r) 3 = 0.47 r] B = 0.55 



(9) 



They all include NLO corrections |) 11]. 

The resulting value for ImA t depends on several input parameters such as \ V c b\ 
\V ub /V cb \ 



m t , B K and JBb^Fb^ Here Bi are the non-perturbative parameters 
related to the hadronic matrix elements of the AS = 2 and AB = 2 operators 
and F Bd is the B d meson decay constant. Values and errors of these input pa- 
rameters, used in the present analysis, are collected in table [5] together with the 
experimental values for e k |35j and Am Bd Except for m t , Er, Am Bd and 



all input quantities are the same as in the review |J7| where further details 



on the chosen ranges with relevant references can be found. See also ||38|| . For 
T{Bd) = 1.6 ps the value for Am Bd in table |3| corresponds to the mixing parameter 
x d = 0.74 ± 0.03 



Table 3: Collection of input parameters. 



Quantity 


Central 


Error 


\Vcb\ 


0.040 


±0.003 


\v ub /v cb \ 


0.080 


±0.020 


B K 


0.75 


±0.15 


\[B~ d F Bd 


200 MeV 


±40 MeV 


MS 

m t 


325 MeV 


±80 MeV 


167 GeV 


±6 GeV 


Am B d 


0.464 ps' 1 


±0.018 ps' 1 


£k 


2.280 • 10" 3 


±0.013- 10" 3 



In what follows we will present two types of the numerical analyses of ImAj 
and e'/e: 

• Method 1: Both the experimentally measured numbers and the theoreti- 
cal input parameters are scanned independently within the errors given in 
table |. 

• Method 2: The experimentally measured numbers and the theoretical input 
parameters are used with Gaussian errors. 
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mean: 1.29 




variance: 0.22 




median: 1.27 




CL: +0.23 


-0.21 (68%) 


+0.49 


-0.38 (95%) 




1.0 1.5 2.0 

\m(\) [1(T 4 ] 



Figure 1: Probability density distribution for lmX t using input parameters as 
given in the text. 



The first method is the one used in our 1993 paper as well as in |J7|. The 
second method is similar to the one used by the Rome group |1| except that 



these authors assume a flat distribution (with a width of 2a) for the theoretical 
quantities. Whereas the first method is even more conservative than adding all 
errors for the various input parameters linearly, the second method yields a similar 
error estimate as if all errors would have been added in quadrature. Thus the 
second method should be considered as reflecting a lower bound on the combined 
error with the true uncertainty lying somewhere in between the two methods. 

In our new analysis let us first concentrate on the case m s (m c ) = 150 ± 
20 MeV. Using the first method and the parameters in table |] we find : 

0.86 • 10~ 4 < ImAi < 1.71 • 10" 4 (10) 

- 1.2 • 10" 4 < e'/s < 16.0 • 10 -4 (11) 

These ranges are similar to the ones found in [[[7] where slightly larger errors for 
m t and Am^ have been used. 

Using next the second method we find the distributions of values for ImA t 
and e' / e in figs. [L] and |2], respectively. From the distributions in figs. [I] and ^] we 
deduce the following results: 

ImAt = (1.29 ±0.22) ■ 10" 4 (12) 
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0.20 



0.15 



CO 

g 



0.05 



mean: 3.6 




variance: 3.4 




median: 3.0 




CL: +3.4 


-2.4 (68%) 


+8.8 


-4.4 (95%) 



0.00 




5.0 10.0 15.0 

e'/e[10 



Figure 2: Probability density distribution for e'/e for m s (m c 
using other input parameters as given in the text. 



150 ± 20 MeV 



e'/e= (3.6 ±3.4) -10" 



(13) 



In addition in the figures we have also given central values according to the median 
prescription and the corresponding 1 and 2 o confidence intervals. Within the 
errors, they agree with mean and variance. We observe that the distribution of 
the values of e' / e is asymmetric with a longer tail towards larger values. However 
at 95% CL. small negative values cannot be excluded. 

The above results for e' / e apply to the NDR scheme, e' / e is generally lower in 
the HV scheme if the same values for B^ 2 ^ and B^ 2 ^ are used in both schemes. 
In view of the fact that the differences between NDR and HV schemes are smaller 
than the uncertainties in B^ 2 ^ and B^ 2 ' we think it is sufficient to present only 
the results in the NDR scheme here. The results in the HV scheme can be found 
in 0,0. 

In spite of some differences in the treatment of hadronic matrix elements 
our results for e'/e, with m s (m c ) = 150 ± 20 MeV, using the second method 
agree well with the results of the Rome group [HJ. On the other hand the 
range in (|TT| ) shows that for particular choices of the input parameters, values 
for e'/e as high as 15 ■ 10~ 4 cannot be excluded at present. Such high values 
are found if simultaneously \V u b/V c b\ = 0.10, B 6 = 1.2, B$ = 0.8, Bk = 0.6, 



UOMeV. 



MS 



405 MeV and low values of m t still consistent with 



Ek and the observed B® — B® mixing are chosen. It is however evident from 
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figure ^] that such high values of e'/e and generally values above 10~ 3 are very 
improbable. 

In [jTS] the hadronic matrix elements relevant for e'/e have been calculated 
within the chiral quark model. Using the first method the authors find a rather 
large range —50 ■ 10 -4 < e'/e < 14 ■ 10~ 4 . In particular they find in contrast 
to |7|, [R| and the present analysis that negative values for e'/e as large as 
—5 • 10~ 3 are possible. This is related to the fact that in the chiral quark model 
Bq and B 8 can for certain model parameters deviate considerably from unity 
and generally B$ > B§. The Dortmund group [jl7|] advocates on the other hand 
B 6 > B 8 to find e'/e = (9.9 ± 4.1) • 10" 4 for m 8 (l GeV) = 175 MeV @. From 



the point of view of the present analysis and the results of the Rome group such 
high values of e'/e for m s (m c ) = 0(150 MeV) are rather improbable within the 
standard model. 



0.20 



mean: 10.4 




variance: 8.3 




median: 8.7 




CL: +9.2 


-5.4 (68%) 


+26.1 


-9.4 (95%) 



>. 

t/5 
c 

<D 
"O 

jg" 0.10 - 
!q 

CO 
-Q 
O 
Q. 



0.05 - 




E'/e [10 ] 

Figure 3: Probability density distribution for e'/e for m s (m c ) = 100 ± 20 MeV 
using other input parameters as given in the text. 



The situation with e'/e in the standard model may however change consider- 
ably if the value for m s is as low as found in [21, 22"fl. Repeating our analysis for 
m s (m c ) = 100 ± 20 MeV we find 



and 



< e'/e < 43.0 • 10~ 4 
e'/e = (10.4 ±8.3) • 10~ 4 



(14) 
(15) 
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in place of (O) and flT3| ) respectively. The corresponding distribution is shown 
in fig. [3[ We observe that the resulting distribution of the values of e' / e is in this 
case rather asymmetric with a very long tail towards substantial positive values. 
Moreover, negative values of e'/e are found to be very unlikely. 

In addition, it is of interest to investigate for which values of the input pa- 
rameters the standard model can reproduce the results of NA31 and E731 collab- 
orations. In table |] we present the values of e'/e for five choices of m s (m c ) and 
for selective sets of other input parameters keeping V cb = 0.040, m t = 167 GeV 
and Bk = 0.75 fixed. The values of e'/e given in this table correspond to sin 5 in 
the first quadrant. The results for the second quadrant turn out to be somewhat 
lower. We observe that the decrease of m s for m s (m c ) > 100 MeV alone is insuffi- 
cient to bring the standard model to agree with the values obtained by the NA31 
collaboration. For central values of other parameters and m s (m c ) = (9(100 MeV) 
the standard model prediction is rather in the ball park of the E731 result. How- 
ever for B 6 > B s , sufficiently large values of | V u b/V c b\ and A^g and small values of 
^s(^c); the values of e'/e in the standard model can be as large as (2 — 4) • 10~ 3 
and consistent with the NA31 result. 

Table 4: Values of e'/e in units of 10~ 4 for specific values of various input 
parameters at m t = 167 GeV, V c b = 0.040 and B K = 0.75. 



\v ub /v cb \ 


A^[MeV] 


B 6 


B 8 


m s (m c ) [MeV] 


e'/e 










75 


16.8 










100 


9.1 


0.08 


325 


1.0 


1.0 


125 


5.3 










150 


3.2 










175 


1.8 










75 


27.8 










100 


15.6 


0.08 


325 


1.2 


0.8 


125 


9.6 










150 


6.2 










175 


4.1 










75 


39.8 










100 


22.5 


0.10 


405 


1.2 


0.8 


125 


14.0 










150 


9.2 










175 


6.2 



To summarize, we have presented a new analysis of e'/e, showing in particular 
the dependence of this important ratio on various input parameters, in partic- 
ular the value of m s . Our results are summarized in (JTT|) and (|i~3"D for the case 
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m s (m c ) = 150±20 MeV and in (0) and QT3[) for the case m s (m c ) = 100±20 MeV. 
Furthermore, table f| gives values of e'/e for particular sets of input parameters. 
The improved results for lm\ t are given in (|T0|) and (|1^) . 

It is clear that the fate of e'/e in the standard model after the improved 
measurement of m t , depends sensitively on the values of \V u b/V c b\, A^jg and in 
particular on B$, B$ and m s . For m s (m c ) = 0(150 MeV) e'/e is generally below 
10~ 3 in agreement with E731 with central values in the ball park of a few 10~ 4 as 
found in [16,37] and here. However, if the low values of m s (m c ) = (9(100 MeV) 



found in [^T],[22| are confirmed by other groups in the future, a conspiration of 
other parameters may give values as large as (2 — 4) ■ 10~ 3 in the ball park of the 
NA31 result. 

Let us hope that the future experimental and theoretical results will be suf- 
ficiently accurate to be able to see whether e'/e ^ and whether the standard 
model agrees with the data. In any case the coming years should be very exciting. 
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